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ABSTRACT 



We have developed an homogeneous evolutionary scenario for H- and He-burning 
low-mass stars by computing updated stellar models for a wide metallicity and age 
range (0.0002< Z <0.004 and 9< t(Gyr) <15, respectively) suitable to study globular 
clusters. This theoretical scenario allows us to provide self-consistent predictions about 
the dependence of selected observational features of the colour-magnitude diagram, 
such as the brightness of the Turn Off (TO), Zero Age Horizontal Branch (ZAHB) 
and Red Giant Branch bump (BUMP), on the cluster metallicity and age. 
f—v ' Taking into account these predictions, we introduce a new observable based on 

^ \ the visual magnitude difference between the TO and the ZAHB [AMy (TO-ZAHB)] , 

(S| ■ and the TO and the RGB-bump [AMy (TO-BUMP)] , given by A = AMy (TO- 

BUMP)-0.566AM y (TO-ZAHB). We show that the parameter A does not depend 
C\J , at all on the cluster age, whereas it does strongly depend on the cluster global metal- 

■ licity. The calibration of the parameter A as a function of Z is then provided, as based 

on our evolutionary models. We tested the reliability of this result by also consid- 

O ering stellar models computed by other authors, employing different input physics. 

1 Eventually, we present clear evidence that the variation of AMy (TO-BUMP) with 

AMy (TO-ZAHB) does supply a powerful probe of the global metal abundance, at 

■ least when homogeneous theoretical frameworks are adopted. Specifically, we show 

that the extensive set of models by VandenBerg et al. (2000) suggests a slightly differ- 
ent calibration of A versus Z calibration, which however provides global metallicities 
higher by only 0.08±0.06 dex with respect to the results from our computations. 

We provide an estimate of the global metallicity of 36 globular clusters in the Milky 
Way, based on our A-Z calibration, and a large observational database of Galactic 
globular clusters. By considering the empirical [Fe/H] scales by both Zinn & West 
(1984) and Carretta & Gratton (1997), we are able to provide an estimate of the 
a— element enhancement for all clusters in our sample. We show that the trend of 
[a/Fe] with respect to the iron content significantly depends on the adopted empirical 
[Fe/H] scale, with the Zinn & West (1984) one suggesting a— element enhancements 
in fine agreement with current spectroscopical measurements. 
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1 INTRODUCTION 

The colour-magnitude diagrams (CMDs) of globular clus- 
ters (GCs) present several features to be compared with the 
various constraints provided by stellar evolution theory. For 
this reason, they are often used to infer fundamental prop- 
erties of ancient stellar populations and, in turn, of the early 
Universe. 

It has been known for a long time that the absolute 
magnitude of the main-sequence turnoff (TO) can be cali- 



brated in terms of the cluster age (t, in Gyr) and chemical 
composition parameters (Y , helium content and Z, global 
metallicity). In order to overcome the uncertainties of ob- 
taining accurate GC distances, it is common to use the 
difference in magnitude between the TO and the zero-age- 
horizontal-branch (ZAHB) at the RR Lyrae instability strip 
(logT e ~3.85), whose luminosity level depends on Y and 
Z, with a negligible dependence on t. Besides this vertical 
method to get the cluster age, there is the so-called horizon- 
tal method, which is based on the determination of the TO 
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colour relative to some fixed point on the red giant branch 
(RGB - see Stetson, Vandenberg & Bolte 1996). As a whole, 
both methods are affected by the empirical difficulties to get 
the precise position of the TO in the CMD, as well as by 
several still-open intrinsic uncertainties on the evolutionary 
models, most notably the dependence on Z of the ZAHB 
luminosity, for the vertical method, and the efficiency of su- 
peradiabatic convection in the outer layers (i.e., the mixing 
length calibration) for the horizontal method, (see Rosen- 
berg et al. 1999, hereafter R99, for details and references). 

The determination of Z for GC stars is by itself prob- 
lematic. Under the classical assumption of scaled-solar metal 
distribution, where all of the abundance ratios of the vari- 
ous metals relative to iron have the solar value, the corre- 
lation between the global metallicity Z and the measured 
iron-to- hydrogen ratio [Fe/H] is given by the classical rela- 
tion logZ=[Fe/H] — 1.70. One should however consider that 
current empirical scales for [Fe/H] (e.g., Zinn & West 1984, 
Zinn 1985, Rutledge, Hesser & Stetson 1997, and Carretta 
& Gratton 1997) may differ by up to ~0.3 dex. In addi- 
tion, it has to be noticed that the assumption of scaled- 
solar chemical compositions might be rather inadequate for 
Galactic GC stars which show a not negligible overabun- 
dance of the so-called a— elements ( O, Ne, Mg, Si, Ar, Ca, 
and Ti) with respect to iron-peak elements ([a/Fe]~ 0.3, 
according to Carney 1996). It follows that the relationship 
between the measured [Fe/H] and the overall metallicity Z 
should follow the relation 

log Z = [Fe/H] - 1.70 + log(0.638/ + 0.362) (1) 

where / is the enhancement factor of a— elements with re- 
spect to iron (Salaris et al. 1993)Q 

The knowledge of the cluster age and global metallic- 
ity is fundamental for studying another important feature of 
CMDs, i.e. the luminosity level of the so-called RGB-bump. 
This is an intriguing feature of the RGB luminosity function 
which has been predicted for a long time, since the pioneer- 
ing theoretical studies by Thomas (1967) and Iben (1968), 
but whose identification in GCs occurred years later (King 
et al. 1985) . A first comparison (Fusi Pecci et al. 1990) of the 
observed difference in visual magnitude between RGB-bump 
and ZAHB level with theoretical models showed that the 
predicted dependence of A1/(ZAHB-BUMP) on metallicity 
was nicely reproduced, but the zero point of the observed 
relation was ~ 0.4 mag too faint. 

An exhaustive investigation on the dependence of 
the parameter AU(ZAHB-BUMP) on the physical inputs 
adopted in stellar computations has been performed by Cas- 
sisi & Salaris (1997, hereinafter CS97). In the same work, it 
was clearly shown that a fine agreement between theory and 
observations does exist when updated RGB stellar models 
and observational data for GCs with spectroscopic measure- 
ments of both iron and a— element abundances, are consid- 
ered. More recently, the same result was obtained by Zoccali 
et al. (1999) and Ferraro et al. (1999, hereinafter F99) by 
using a much larger sample of Galactic GCs than CS97. 

* It is worth recalling that this relation holds strictly only in 
case of the scaled-solar metal mixture by Ross & Aller (1976). 
The coefficients of this relation has to be slightly changed when 
accounting for different solar metal distributions. However, these 
changes are very small (see Yi et al. 2001). 



However, one has also to notice that Bergbusch & Vanden- 
berg (2001) have recently found a discrepancy of the order 
of 0.25 mag, when comparing the AV(ZAHB-BUMP) values 
provided by their models with the observational data for 4 
GCs at various metallicities. 

In reality, the RGB-bump luminosity is also dependent 
on the age of the stellar population (<5V(BUMP)/<51ogt ~ 
1), while the above investigations do not take into ac- 
count the age of individual GCs. It is well known the cur- 
rent debate whether the Galactic GCs formed all at the 
same time or there was a significantly protracted forma- 
tion epoch. A recent analysis of an homogeneous photomet- 
ric database containing 34 Galactic GCs (R99) shows that 
within the intermediate-metallicity clusters ( — 1.2 <[Fe/H]< 
—0.9) there is a clear evidence of age dispersion, with some 
clusters up to ~30% younger than the oldest ones. As an ex- 
ample, we consider the globular cluster NGC 362 for which 
F99 give AF(ZAHB-BUMP)=0.10 mag±0.12 mag and a 

global metallicity0 [M/H]=-0.99. 

From the relations given by these authors 

M v ,bump = 0.75+0.99 log t+1.58[M/H]+0.26[M/H] 2 (2) 

My, ZAHB = 1.00 + 0.35[M/if] + Q.05[M/H] 2 (3) 

one has that the agreement between theory and observa- 
tions occurs for a cluster age of ~14.9 Gyr, which is the 
value derived by R99 for the bulk of coeval GCs, using the 
Straniero et al. (1997) theoretical isochrones. However, R99 
find that the actual age of NGC 362 is ~ 11.5 Gyr, in which 
case one obtains a AV(ZAHB-BUMP) value by 0.21 mag, 
namely ~0.1 mag larger than the observed one. 

It should be noticed that R99 adopts the [Fe/H] scale 
by Carretta & Gratton (1997), without any correction due 
to a— enhancement and, in order to interpret the observed 
difference AU(TO-ZAHB), a linear relation for the abso- 
lute magnitude of the ZAHB as a function of the metal 
content (Mv,zAHB=0.18([Fe/H]+1.5)+0.65), instead of the 
quadratic relations used by F99 [see Eq. (3)] and CS97 
(M v ,zahs = 1.13+0.39[M/H]+0.06[M/H] 2 ). In other words, 
the quoted extensive investigations on the RGB-bump lumi- 
nosity and the age of Galactic GCs are not mutually con- 
sistent, since they use different metallicity scales and adopt, 
as theoretical values for the RGB-bump and TO luminosity, 
different calibrations of the ZAHB luminosity in terms of 
metallicity. 

In order to make a complete and self-consistent study of 
these most prominent features in CMDs, we present in sec- 
tion 2 an homogeneous set of theoretical evolutionary results 
which predict the TO, RGB- bump and ZAHB luminosity for 
a wide range of ages (9< t(Gyr) <15) and metal abundances 
(0.0002< Z <0.004). In section 3 we introduce a new param- 
eter A defined as A = AMv(TO-BUMP)-0.566AM v (TO- 
ZAHB), and we show that this depends only on the global 
metallicity Z, being unaffected by age. In order to test 
this result, the evolutionary models computed by Salaris 
& Weiss (1997, 1998) and by Vandenberg et al. (2000) 
are also considered. We show that the former ones are in 
close agreement with our computations, whereas in the lat- 
ter case some significant differences are present, leading to 

t [M/H] = [Fe/H]+log(0.638/+0.362). 
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Table 1. Selected parameters from the evolutionary models. 
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a slightly different formulation of the AMy (TO-BUMP) 
versus AMy (TO-ZAHB) relationship. Nevertheless, the pa- 
rameter A—AMv (TO-BUMP) —0.487AMy (TO-ZAHB) in- 
ferred from the Vandenberg et al.'s (2000) computations is 
again strongly dependent on global metal abundance with 
no dependence on age. 

On this basis, the TO, ZAHB and RGB-bump observa- 
tions collected by F99 and R99 for a large sample of Galactic 
GCs are interpreted. The resulting global metal abundances 
are presented in section 4, together with some clues to the 
a— element enhancement with respect to iron in Galactic 
GCs. A summary is finally given in the last section. 



2 THEORETICAL MODELS 

In order to understand the observed properties of CMDs and 
RGB luminosity functions of GCs, we need an homogeneous 
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Figure 1. Theoretical predictions for AMy (TO-BUMP) and 
AMy (TO-ZAHB) as a function of age (t) and global metallic- 
ity (Z). The dashed lines refer to constant ages in the range of 
9 to 15 Gyr (left to right), while the solid lines refer to constant 
mctallicitics in the range of 0.0002 to 0.004 (bottommost line). 
The solid lines have a similar slope of 0.566 (see text). 

set of evolutionary models for both the H- and He-burning 
phases. In present work, we have adopted a set of theoretical 
models computed by using the FRANEC evolutionary code 
(Chieffi & Straniero 1989; CS97; Castellani et al. 1997). A 
subset of these stellar models and corresponding isochrones 
have been already presented by Cassisi et al. (1998, 1999). 
However, in order to improve the coverage of the wide metal- 
licity range spanned by Galactic GCs, we have extended 
previous computations considering additional initial chemi- 
cal compositions. The physical inputs adopted in the model 
computations have been already discussed by Cassisi et al. 
(1998). A short summary is given below. 

We use the OPAL (Rogers et al. 1996) equation of state 
(EOS), supplemented, in the temperature-density regime 
not covered by the OPAL tables, by the EOS provided by 
Straniero (1988), plus a Saha EOS in the outer stellar layers. 
The radiative opacities provided by the OPAL group (Igle- 
sias, Rogers & Wilson 1992) for temperatures larger than 
10,000K are employed, together with the Alexander & Fer- 
guson (1994) tables for lower temperatures. Electron con- 
duction opacities come from Itoh et al. (1983). The bound- 
aries of the convective regions are determined according to 
the canonical Schwarzschild criterion. Induced overshooting 
and semiconvection during the central He-burning phase 
have been accounted following the prescriptions given by 
Castellani et al. (1985). The standard mixing- length theory 
is used for evaluating the temperature gradient in the su- 
peradiabatic regions; the value of the mixing length param- 
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logZ 

Figure 2. (top) - Mean value of the parameter A=AMv(TO- 
BUMP)-0.566AM V (TO-ZAHB) as a function of Z. The solid 
line is the analytic fit to the data (see Eq. (4) and Eq. (5) in the 
text), (bottom) - Difference between the evolutionary metallicity 
logZ ell derived using Eq. (4) and Eq. (5) and the original value 
logZ, as a function of logZ. The discrepancy is within ±0.03 dex 
(dashed line). 



eter is fixed at 1.6 for all computations. This value allows 
to match the RGB effective temperatures of a selected sam- 
ple of Galactic GCs provided by Frogel, Persson & Cohen 
(1983). 

All the models include microscopic diffusion of both he- 
lium and heavy elements, using the formalism by Thoul et 
al. (1994). For a detailed discussion concerning the effects 
of atomic diffusion on both main-sequence stars and more 
evolved models, we refer the interested reader to Castel- 
lani et al. (1997) and Cassisi et al. (1998). ft is well known 
that helioseismology provides clear evidence that helium and 
heavy element settling is occurring in the Sun (see Basu, 
Pinsonneault, & Bahcall 2000, and references therein). How- 
ever, recent spectroscopic analyses of GC stars (Gratton et 
al. 2001; Castilho et al. 2000; Ramirez et al. 2001) show 
that the surface chemical patterns do not show the effects of 
atomic diffusion, at least for what heavy elements are con- 
cerned. This means that either chemical settling is less effi- 
cient in metal-poor stars, or some other physical mechanism 
at work in the stellar atmospheres counterbalances the ef- 
fect of atomic diffusion. A firmer understanding of this topic 
has to wait for more accurate theoretical and observational 
investigations. 

Evolutionary models with masses ranging from 0.6 to 
1.3M© have been computed from the Zero Age Main Se- 
quence up to the tip of the RGB marking the triple-alpha 
ignition inside the degenerate He core. The explored range of 



metallicity is 0.0002 < Z < 0.004, while for each metallicity 
the same initial He abundance F=0.23 has been adopted. 
The He-core mass M c and the chemical stratification at the 
RGB He-flash provided by the different evolutionary models 
have been used to compute the ZAHB structures needed to 
estimate the ZAHB luminosity level at logT' e =3.85, the aver- 
age effective temperature of the RR Lyrae instability strip, 
fsochrones have been then constructed in the age range 9 - 
15 Gyr. 

Finally, all the models and isochrones have been trans- 
posed from the theoretical H-R plane to the observational 
(V, B — V) plane by using the bolometric corrections and 
colour-temperature relations provided by Castelli, Gratton 
& Kurucz (1997a,b). In Table 1, for the different chemi- 
cal compositions adopted in present work, we report the 
TO, RGB-bump and ZAHB absolute visual magnitude, as 
a function of age. 



3 THEORETICAL CALIBRATIONS 

The data listed in Table 1 show that, at a fixed age, TO, 
RGB-bump and ZAHB luminosities decrease with increasing 
the metal content. At a fixed metallicity the TO and RGB- 
bump become fainter for larger ages, whereas the ZAHB 
level becomes slightly brighter, due to the slight increase 
of the He-core mass M c with decreasing total mass of the 
progenitor on the RGB. 

Figure 1 displays the two differential parameters 
AMv (TO-ZAHB) and AMy (TO-BUMP), as derived from 
the results in Table 1. The dashed lines show the theoretical 
trend at a constant age, in steps of 1 Gyr, from 9 Gyr (left) 
to 15 Gyr (right), while the solid lines refer to a constant 
Z from 0.0002 to 0.004 (bottommost line). The difference 
AMv(TO-BUMP) increases for larger ages, at fixed Z, and 
decreases for larger metallicities, at fixed t. As for the pa- 
rameter AMv (TO-ZAHB), it increases with increasing Z up 
to logZ ~ —2.7, and then it decreases for larger metallici- 
ties, at least when Z <0.004. On this subject, we recall that 
all our models are computed with an initial helium abun- 
dance y=0.23. Within the quoted range of metal content, 
the expected He increase due to an helium-to-heavy element 
enrichment ratio AY/AZ ~2.5 is only AY=0.01, whereas 
for metallicities larger than Z=0.004 the original amount of 
helium is expected to become significantly larger. Further- 
more, also the amount of extra-helium brought to the stellar 
surface by the first dredge up is expected to increase when 
moving towards larger metallicities. As a consequence, the 
above theoretical trends cannot be applied to very metal-rich 
clusters ([Fe/H]> -0.6). 

The data in Fig. 1 disclose that, for ages between 9 and 
15 Gyr, the variation of AMy (TO-BUMP) with AMv (TO- 
ZAHB) has a constant slope of 0.566 (solid lines), with 
the zero-point depending on Z. From the results listed in 
Table 1 we calculate then the mean value of the parame- 
ter A— AMv (TO-BUMP) — 0.566AMy (TO-ZAHB) for each 
given Z. The upper diagram in Fig. 2 shows that the results 
are well approximated by the following linear relations: 

log Z = +0.518 - 1.933A (4) 
© 0000 RAS, MNRAS 000, 000-000 
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Figure 3. (top) - Comparison between our theoretical predictions Figure 4. As in Fig. 3, but for the TO luminosity, 

for the RGB Bump and ZAHB luminosity and the ones provided 
by Salaris & Weiss (1997, 1998: open and filled circles, respec- 
tively) and Yi et al. (2001: asterisks). Each point corresponds to 
a fixed age and metallicity, for the values in common with our 
models, bottom) - As top panel, but for the models by Vanden- 
berg et al. (2000). Open circles refer to the ZAHB luminosity, 
while full circles correspond to the RGB Bump luminosity. 



logZ = -1.110 - 1.039,4 



(5) 



with A > 1.82 and A < 1.82, respectively. As it is shown in 
the lower panel of Fig. 2, the discrepancy between the evo- 
lutionary estimate of the global metallicity (hereafter Z ev ) 
based on the two differential parameters AMy (TO-BUMP) 
and AMy-(TO-ZAHB), and the original value, is close to 
±0.03 dex. 

Before proceeding further, it is worthwhile to com- 
pare our evolutionary predictions with similar computa- 
tions provided by different authors. We considered the re- 
cent isochrones and ZAHB models computed by Salaris & 
Weiss (1997, 1998, hereinafter SW) and by Vandenberg et 
al. (2000, hereinafter Vdb, see also Bergbusch & Vandenberg 
2001) with Z ranging from 0.0002 to 0.004 and ages in the 
range of 10-14 Gyr and 10-16 Gyr, respectively. We wish to 
notice that both sets of models do not account for atomic 
diffusion and have been computed by adopting a slightly 
different AY/AZ ratio (3 and 2, respectively). Moreover, 
they are computed under different assumptions about the 
a— element enhancement. In our comparison we used the 
Vdb models based on a scaled-solar heavy element distribu- 
tion. 

As shown in the upper panel of Fig. 3, both the ZAHB 
(open circles) and RGB-bump (filled circles) absolute mag- 
nitudes inferred from SW models are in close agreement 
with our results, whereas the Vdb computations plotted 
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Figure 5. As in Fig. 1, but in this case the Vdb models have 
been adopted. 
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Figure 6. As in Fig. 2, but adopting the Vdb models. 



in the lower panel suggest brighter RGB-bump and fainter 
ZAHB luminosities, with respect to the data in Table 1. The 
upper panel of Fig. 3 shows also the comparison between 
the RGB-bump absolute magnitudes provided by the recent 
Yale isochrones (Yi et al. 2001) and our results (asterisks), 
for the common values of metallicity (Z=0.001 and 0.004) 
and age (10 and 13 Gyrs). As a whole, one has that the Yale 
isochrones suggest brighter RGB-bumps than our models. 
Note that we were unable to carry out a similar comparison 
with the Yale ZAHB absolute visual magnitudes, as no data 
have been provided in the literature. Figure 4 deals with 
the TO luminosity provided by the quoted sets of computa- 
tions, compared with the results in Table 1. We get that SW 
models are again in agreement with our predictions, whereas 
VdB and Yale models give brighter absolute magnitudes. 

On the basis of the above discussion, we have considered 
the models by VdB for testing the effect of different input 
physics on the scenario suggested by our computations. Fig- 
ure 5 presents AMy-(TO-BUMP) versus AM^ (TO- ZAHB), 
as derived from Vdb models. The dashed lines show the the- 
oretical trend at a constant age, in steps of 2 Gyr, from 10 
Gyr (left) to 16 Gyr (right), while the solid lines refer to 
a constant Z from 0.0002 to 0.004 (bottommost line). The 
general trend in Fig. 1 is here confirmed, but the variation 
of AMy(TO-BUMP) with AMv(TO-ZAHB) has a slope of 
0.487 (solid lines). The correlation between the mean value 
of the parameter A=AMy (TO-BUMP)-0.487AAfv (TO- 
ZAHB) and logZ is presented in the upper diagram of Fig. 
6, where the solid line depicts the relations 

log Z = +1.617- 2.119,4 (6) 



Figure 7. {top) - The evolutionary global metallicity logZ e „ for 
the GCs in Table 3, as a function of the measured [Fe/H] value 
(with an uncertainty of ± 0.10 dex) on the Carretta & Gratton 
(1997) scale. The solid line represents the canonical correlation 
for scaled-solar chemical compositions [/=1 in Eq. (1)], while the 
dashed line refers to /=3. (bottom) - The a— enhancement fac- 
tor / for the clusters in Table 3 as a function of the measured 
[Fc/H] value on the Carretta & Gratton (1997) scale. The solid 
line represents the classical value for scaled-solar chemical com- 
positions (/=1), while the dashed lines denote the average values 
/ ~ 1.5 and / ~ 3 suggested by clusters with [Fe/H]> —1.7 and 
[Fe/H]< -1.7, respectively. 

log Z= -0.154 - 1.330A (7) 

with A > 2.24 and A < 2.24, respectively. As shown in 
the lower panel of Fig. 6, also in this case the discrepancy 
between the evolutionary estimate of the global metallicity 
(hereafter Z ev ) and the original value is close to ±0.03 dex. 



4 GLOBAL METALLICITY FOR GLOBULAR 
CLUSTERS IN THE GALAXY 

In this section we provide accurate evolutionary estimates 
of the global metallicity of a large sample of Galactic GCs, 
by using the relationships discussed in the previous section. 

The sample of Galactic GCs studied in this paper is 
listed in Table 2. The [Fe/H] values are listed both for the 
Zinn & West (1984) and the Carretta & Gratton (1997) 
metallicity scale, while the TO visual magnitudes are from 
R99, except for the following clusters: NGC 4833 and 5286 
(Samus et al. 1995a,b), NGC 5694 (Ortolani & Gratton 
1990), NGC 6584 (Sarajedini & Forrester 1995), NGC 6717 
(Ortolani et al. 1999), NGC 6934 (Piotto et al. 1999) NGC 
7006 (Buonanno et al. 1991) and NGC 7492 (Buonanno et 
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Table 2. Observed parameters for Galactic Globular Clusters. Iron-to-hydrogen contents are given in the Zinn & West (1984, [Fe/H]#w) 
and Carretta & Gratton (1997, [Fe/HJ^g) empirical scales. For the latter scale, the fourth column gives the source: = Carretta & 
Gratton (1997); 1 = Rutledge, Hesser & Stetson (1997); 2 = Ferraro et al. (2000); 3 = Rosenberg et al. (1999). The [Fe/H] CG value of 
NGC 7492 has been estimated through the mean relation between [Fe/HJ^c and [Fe/H]^;,^. 
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al. 1987). The ZAHB and RGB- bump observations are from 
F99, except for the following additional measurements of the 
RGB-bump level: NGC 6352 and 6397 (Alves & Sarajedini 
1999) and NGC 6362 (Piotto et al. 1999). We wish to no- 
tice that the two major sets of observations are mutually 
consistent as the ZAHB visual magnitudes listed by F99 are 
~ 0.07 mag fainter than the HB visual magnitudes in R99, 
except for NGC 6093 and NGC 6254. For these two clus- 
ters, the F99 values are brighter than those given by R99 
(see later). Moreover, we have also checked that the ZAHB 
magnitudes given by F99 are consistent with the additional 
measurements of the TO and RGB-bump visual magnitudes. 

As a first step, we determine the values of A=AV(TO- 
BUMP)-0.566A1/(TO-ZAHB) from the data in Table 2, 
then we derive the global metallicity using Eq. (4) and Eq. 
(5). The results for Z ev are listed in column (4) of Table 
3, together with the horizontal branch morphological type 
(HB-type) | from Harris (1996). For the two clusters NGC 



t The HB-type is given by the ratio (B-R)/(B+V+R), where B 
and R denote the number of HB stars blueward and redward of 



6093 and 6254, the last two lines in Table 3 refer to the HB 
level given by R99. 

The correlation between \ogZ ev and the measured 
[Fe/H] parameter^ is presented in the upper panel of Fig. 
7, where the solid line shows the canonical correlation 
logi?= [Fe/H]— 1.7 for scaled-solar chemical compositions. 
One obtains that almost all the clusters present a mild over- 
abundance of the global metallicity with respect to the mea- 
sured iron-to- hydrogen content. In terms of a— element en- 
hancement [see Eq. (1)], we show in the bottom panel in Fig. 
7 that the clusters with [Fc/H]> —1.7 show / ~ 1.5, whereas 
the four most metal-poor clusters show / ~ 3 (dashed lines) . 

Figure 8 displays the results based on Vdb computa- 
tions. These models suggest slightly larger Z av values (see 
data listed in column (5) of Table 3) in comparison with our 
models, with an average difference of 0.08±0.06 dex which 
is within the uncertainty of the derived global metallicity. 

the RR Lyrae instability strip, respectively, and V is the number 
of RR Lyrae stars. 

§ Note that the measured [Fe/H] values are taken with an uncer- 
tainty of ±0.10 dex. 
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-3.5 




-2 -1.5 
[Fe/H] cl 




Figure 8. As in Fig. 7, but for Vdb models. The dashed line in 
the top panel refers to /=3.7, while the dashed lines in the bottom 
panel denote the average values / ~ 1.9 and / ~ 3.7 suggested 
by clusters with [Fe/H]> —1.7 and [Fe/H]< —1.7, respectively. 




Figure 9. As in Fig. 7, but for the [Fc/H] scale by Zinn & West 
(1984). The dashed line in both panels refer to an a— enhancement 
factor / ~ 3 suggested by clusters with [Fc/H]< —1.3. 



Figure 10. As in Fig. 9, but for Vdb models. The dashed line in 
both panels refers to an a— enhancement factor / ~ 3.5 suggested 
by clusters with [Fc/H]< -1.3 



Also the a— element enhancement is slightly increased, clus- 
ters with [Fe/H]> —1.7 now showing / ~ 1.9, whereas the 
four most metal-poor clusters show / ~ 3.7 (dashed lines). 

We remark here that the global metallicity given in Ta- 
ble 3 are derived from the differential parameters AV(TO- 
ZAHB) and AV(TO-BUMP) and are fully independent of 
the measured iron-to-hydrogen content. On the contrary, the 
/ enhancement factors plotted in the lower panel of Figs. 7 
and 8 obviously depend on the adopted empirical [Fe/H] 
scale. When we consider the Zinn & West (1984) scale the 
results are quite different. As shown in Fig. 9 and Fig. 10, 
the a— enhancement factor is / ~ 3 (our models) or ~ 3.5 
(Vdb models) with [Fe/H]< —1.3, declining toward the solar 
value / ~1 at larger metallicities. 

The abundances of a— elements in Population II stars 
provide fundamental constraints to the chemical evolution 
models of the Galaxy. Several spectroscopic studies have 
been devoted to this problem (see Carretta, Gratton & 
Sneden 2000 for discussion and references). As a whole, 
a— elements appear to be overabundant by a constant factor 
~2-3 in metal-poor stars, declining toward scaled solar ra- 
tios around [Fe/H]=-1.0. As the run of [a/Fefl with [Fe/H] 
is required to study the chronology of the Galactic halo for- 
mation (Matteucci & Francois 1992), we notice that mild 
variations of the empirical [Fe/H] scale may lead to quite 
different pictures for the dependence of [a/Fe] on [Fe/H], 
and different [Fe/H] values for the transition toward scaled- 
solar ratios. 



[a/Fe]=log/. 
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Table 3. Global metallicity for our selected sample of Galactic 
GCs. 



NGC 


[Fe/HW 
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Both the distribution and absolute values of the 
a— element enhancement suggested by our method together 
with the Zinn & West (1984) [Fe/H] scale appear in fine 
agreement with the empirical results by Carney (1996). On 
the contrary, the results obtained using the Carretta & Grat- 
ton (1997) [Fe/H] scale do not appear fully supported by 
current spectroscopic measurements. Even though to inves- 
tigate the accuracy of current empirical [Fe/H] scales is out 
of the aim of the present paper, we wish to emphasize that 
this is still a controversial issue. We also wish to notice 
that recently Vandenberg (2000) and Bergbusch & Vanden- 
berg (2001) have provided some additional evidence, based 
on the comparison between their theoretical isochrones and 
CMDs of Galactic GCs, that the metallicities determined 
by Carretta & Gratton (1997), mainly for the intermediate- 
metallicity clusters, are too high, values close to the Zinn 
& West (1984) scale being more favourite. Present analysis 
provides independent support to their claim. 

Before closing this section, it is worthwhile to emphasize 
that our method, being based only on well defined observa- 
tional features of the CMD, appears quite suitable for deriv- 
ing accurate estimates of the GC global metallicity. It uses 
only differential quantities, i.e. visual magnitude differences, 



and its results are therefore insensitive to uncertainties on 
the zero point of the CMD photometric calibration as well 
as on the cluster distance and reddening. 

It is also worth noticing that the method relies on ob- 
servables whose properties are not significantly affected by 
non-canonical processes occurring along the RGB, like deep- 
mixing processes, which are currently invoked to explain the 
chemical abundance patterns observed in RGB stars of many 
Galactic GCs (Sweigart & Mengel 1979, Langer, Hoffman & 
Sneden 1993). As suggested by theoretical predictions (Mes- 
tel 1957) and spectroscopic data (Charbonnel, Brown, & 
Wallerstein 1998, and references therein), molecular weight 
gradients, such as the H-discontinuity left over in the en- 
velope of RGB stars by the first dredge up, should be able 
to inhibit these non-canonical mixing processes, which are 
therefore expected to be efficient only after the RGB bump. 



5 SUMMARY 

The main results obtained in present work can be summa- 
rized as follows: 

• We have developed an homogeneous theoretical sce- 
nario, based on updated stellar models including atomic dif- 
fusion, computed for a wide metallicity range; this allows us 
to predict the behaviour with the cluster metallicity and age 
of relevant observational features of the CM diagram, such 
as the Turn Off magnitude and the brightness of both the 
ZAHB at the level of the RR Lyrae instability strip, and the 
RGB-bump; 

• we define a new observable based on the visual magni- 
tude difference between the TO and the ZAHB AMy(TO- 
ZAHB), and the TO and the RGB-bump AM v (TO- 
BUMP), namely, A = AM v (TO-BUMP)-0.566AMy (TO- 
ZAHB). We show that the parameter A does not depend 
at all on the cluster age, whereas it does strongly correlate 
with the global cluster metallicity. The calibration of the pa- 
rameter A as a function of Z is provided; 

• by using this calibration, together with a large obser- 
vational database of Galactic GCs, we estimate the global 
metallicity of a large sample of GCs. Our estimates are inde- 
pendent of empirical metallicity measurements, but strongly 
rely on the accuracy and reliability of the adopted theo- 
retical models. By using our global metallicity determina- 
tions, together with presently available [Fe/H] determina- 
tions (Zinn & West 1984 or Carretta & Gratton 1997), and 
using the relation provided by Salaris et al. (1993) that con- 
nects global metallicity to [Fe/H] and a— element enhance- 
ment, we are able to provide an estimate of the [a/Fe] ratio 
for the individual clusters in our sample; 

• the trend of [a/Fe] with [Fe/H] we predict when using 
the Zinn & West metallicity scale, appears in fine agreement 
with current spectroscopical measurements, which suggest 
a rather constant overabundance of a— elements in metal- 
poor stars, and a decline toward the solar ratios around 
[Fe/H]~ —1.0. This picture is not supported when we use 
the Carretta & Gratton [Fe/H] scale; 

• the calibration of the A — Z ev relation is model de- 
pendent, in the sense that the Vdb models, which are 
based on different input physics, yields a slightly different 
formulation of the parameter A, given by A=AM^(TO- 
BUMP)-0.487AMv(TO-ZAHB). However, the resulting 
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Ze V values are only 0.08±0.06 dex larger than the results 
based on our computations; 

• the good agreement between the global metallicity esti- 
mates provided by our models accounting for atomic diffu- 
sion, and the ones obtained by employing Vdb models which 
neglect atomic diffusion, shows that the theoretical calibra- 
tion of the A parameter is only marginally affected by the 
efficiency of this non-canonical process; 

• the present analysis clearly demonstrates that the our 
newly defined observational parameter, independent of the 
uncertainty on the zero point of the photometric calibration 
as well as on the cluster distance and reddening, is a powerful 
indicator of globular cluster global metal abundances. 
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